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Abstract 

The optical properties of wide Quantum Wells are considered, taking into account the screened 
electron-hole interaction potential and parabolic confinement potentials, different for the electrons and 
for the holes. The role of the interaction potential which mixes the energy states according to different 
quantum numbers is stressed. The results obtained by our method are in agreement with the observed 
spectra and give the possibility to the assessment of the resonances. 


1 Introduction 

A single Quantum Well (QW) is formed when a thin layer of a narrow gap semiconductor material (Well) 
lies between layers of wider gap materials (Barriers). The size of QWs (taken in the growth direction) is 
generally between 1 and 20 nm (narrow QWs). Mostly the size is comparable with the effective excitonic 
Bohr radius a* of the QW material (for example, a* is about 10 nm for GaAs). When the size of the QW is of 
the order of a few excitonic Bohr radii, we speak about Wide Quantum Well (50-100 nm in the case of GaAs 
Wells). Typical property of the QWs is the confinement of the carriers (electrons and holes) inside the well. 
Mostly one considers a rectangular shape of the confinement potential, but also other shapes are considered. 
A special fabrication technique allows to create a parabolic shaped confinement, leading to the so-called 
Parabolic Quantum Wells (PQWs). More generally, we can speak about Parabolic Quantum Nanostructures 
(for example. Parabolic Quantum Dots). When the size of a PQW corresponds to the Wide QW region, 
we deal with Wide Parabolic Quantum Wells (WPQWs). The structures with parabolic confinement have 
attracted more attention in the recent decades (for example. Ref. [I]-[T5]'). 

Each type of semiconductor structures has a specific property depending on its dimensionality, which 
influences the optical characteristic of a given structure. In bulk crystals the coupling between the external 
electromagnetic wave and the internal polariton modes gives rise to the so-called ABG problem, which was 
extensively examined in the past [l6]-[25]). In low-dimensional structures as Quantum Wells, Quantum Wires 
and Quantum Dots mostly the long-wave approximation is used, so the electrodynamical aspect is simplihed, 
however, on the other hand, due to the confinement, the separation of center-of-mass and relative e-h motion 
is not possible and the problem of solving 6-dimensional Schrodinger equation appears. The additional 
difficulty consists on the fact that the spherical symmetry of the Coulomb potential is not compatible with 
the cylindrical symmetry of most nanostructures. Thus several approximations have been proposed which 
used the lowering of the dimension of the problem (for example, [IS]). Recently, also due to the increasing 
power of computers, direct numerical solution of the 6-dimensional 2-particles Schrodinger equation has been 
performed (for example, Refs. ^7}, [28]). 

Here we consider Wide Quantum Wells (WQWs), nanostructures where the optical active layers have 
extension of a few excitonic Bohr radii in the growth direction. To a certain approximation, the polaritonic 
aspect can be here neglected, and the long-wave approximation sustained. The exciton is not squeezed as 
in other low-dimensional structures. In typical Quantum Wells with the dimension of, say, one excitonic 
Bohr radius in the growth direction we observe only a few excited states. In WQWs, due to the greater 
extension, significantly larger number of states is observed. The Coulomb potential and different confine¬ 
ment potentials for electrons and holes couples electron and hole confinement states of different quantum 
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numbers. Such phenomena have been observed experimentally (see, for example, Ref. m)- We propose the 
computational method which leads to analytical expression for the electric susceptibility of wide parabolic 
quantum well taking into account the screened electron-hole interaction and parabolic confinement potential. 
With the purpose of exemplification, we consider a quantum well with GaAs as the optically active layer 
and Gai_a;Ala;As as the barriers, where the active layer is of the extension of a few excitonic Bohr radii. 
The absorption spectra of such a structure show a large number of resonances (n = 8 observed in my The 
choice of optimal effective potential parameters as well as the damping constant used in our calculation is 
verified by numerical calculations of the total fitting error for maxima of susceptibility. We have chosen 
as reference the paper by Miller et al [T] because it contains a lot of experimental data which allowed to 
compare the obtained theoretical results with experiment. The agreement between our calculated spectrum 
and experimental data is very good with regard to the number and position of the maxima of susceptibility. 

Our paper is organized as follows. In the section 2, we present the assumptions of considered model 
and solve the constitutive equation with effective electron-hole interaction potential. Section 3 is devoted to 
the details of the applied potential. Next, in section 4, the derived solution of constitutive equation is used 
to obtain the energy levels of the considered GaAs/Gai_a;Al 3 ,As wide parabolic quantum well. Finally, in 
section 5, the susceptibility for such nanostructure is calculated and discussed. The comparison of obtained 
results with experimental data and a brief overview of optimizing procedure is included. 


2 The Model 


We will compute the linear optical response of a WPQW to a plain electromagnetic wave 

Ei{z,t) = ifioexp(ifco2: - iwf), fco = —, (1) 

c 

attaining the boundary surface of the WPQQW active layer located at the plane z = 0. The second boundary 
is located at the plane z = L. In the case of GaAs Well the extension L will be of the order 2-40 nm. 

Due to this extension, the following aspects should be taken into account. I) Several confinement states 
resulting from the confinement in the z direction are to be included in the consideration. 2) The electron- 
hole potential plays an important role and cannot be approximated by a 2-dimensional potential, as was 
sometimes done in the case of simple Quantum Wells (the limit L —>■ 0). 3) The parabolic shape of the 
confinement potential is assumed 


Vu^iZe,Zh) = -b ^mhzColzl. 


( 2 ) 


4) The electron-hole interaction is described by the potential V (re, r/j). 5) We adopt the real density matrix 
approach to compute the optical properties. In this approach the linear optical response will be described by 
a set of coupled equations: two constitutive equations for the coherent amplitudes W(re, r?i), = II,L stands 

for heavy-hole (H) and light-hole exciton); from them the polarization can be obtained and used in Maxwell’s 
field equations. Having the field we can determine the QW optical functions (reflectivity, transmission, and 
absorption). 

Thus the next steps are the following: We formulate the constitutive equations. The equations will be 
then solved giving the coherent amplitudes Y. From the amplitudes we compute the polarization inside the 
Quantum Well, the electric field of the wave, and the optical functions. This scheme will be applied for the 
case investigated in Ref. [T]. 

As was explained in, for example. Ref. |24j . the constitutive equation for the coherent amplitude T in a 
Quantum Well has the form 


En — hu) — ir + 


pIz 


+ 


pIz 


Pp 


“b /f ^eh Zq^ Zfi) -b Yconi(^Z q^ Z/j) 


2me ‘2mhz 2/X|| 2M|| 


Y = M(r)E(R) 


(3) 


where M(r) is the transition dipole density, which form we have assumed as 

M(r) = M(p, z, 4 ,) = (p _ pg) ^ 

27rpo 


( 4 ) 
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z = Ze — Zh being the relative coordinate in the z direction, po is the coherence radius, R jest is the 
excitonic center-of-mass coordinate and E(R) is the electric field vector of the wave propagating in the QW; 
Vconfize, Zh) is the confinement potential for electrons and holes, and Pp, p|| are the momentum operators 
for the excitonic relative- and center-of-mass motion in the QW plane. 

In the following we assume that the propagating wave is linearly polarized in the x direction, and that 
the vector M has a non-vanishing component in the same direction. Taking the confinement potential in the 
form ([2]) we find in the equation ((S]) Hamilton operators for the one-dimensional harmonic oscillator 


R. = 


2m' ^ 2 ”"'="''^'’ 


Hh = 


pIz 

2 mh 


-mhUjlzl. 


(5) 


Therefore we look for a solutions Y in terms of the eigenfunctions of the operators He^Hh 


N 

Y{p,Ze,Zh)= 'tpej{Ze)'tl^nh{zh)Yjn{p), ( 6 ) 

j,n=0 


The eigenfunctions tpj have the form 


IpejiZe) = TT 


= ^-1/4 


21 j\ 


-Hj{aeZe)e 2 "^=; ae = 


niezUJe 


mhzUJh 


fpnhizh) = tt ^^Hn{ahZh)e 2 '“^; ah = 

with the Hermite polynomials Hn{x), and the corresponding eigenvalues + 5 ) Substituting 

into the eq. ([3]) we obtain equations for the functions Yj^ 


(7) 
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E 

j,n=0 


Pp Py 

Eg — hu! — ir + Eje + Enh + 7 ^- 1 - „ n,f + Yeh {P, Ze, Zh) 

2p|| 2M|| 


ipjize)i^n{zh)Yjn{p) = M(r)E(R). 


( 8 ) 


Now we have to specify the shape of the interaction potential Veh{p, Ze, Zh) and the wave electric field E(R). 
We assume the so-called long-wave approximation and consider E(R) in the equation (|S]) as a constant 
quantity.The electron-hole interaction potential Veh{p, Ze, Zh) is, in general, the screened Coulomb potential 


Veh{p, Ze, Zh) — , -, (9) 

47re&Yp2 _|_ _ 2.^)2 

£{, being the dielectric constant of the QW material. Despite of the nanostructures with cylindrical symmetry 
considered in ref. [28], in the case of the wide QWs one does not have an orthonormal basis of functions so 
the use of an effective e-h interaction potential will be made 


Veh 


— S exp 




( 10 ) 


where r;, w are certain parameters which will be estimated below. Using the above potential, the dipole 
density (jd]), and neglecting the center-of-mass in plane motion, we put the constitutive equation ([3|) into the 
form 


Er 


) b;, - V VrsnjYr^g = E^{r\s)5 {p - Po) ; 

2pi\ J ^ 2ttpo 


( 11 ) 


where 


Ers = Eg + Ere + Esh - - lE, r, s, = 0,1,2,... 


Vr 


rsnj 


= S{r 


exp 


-V {Ze - Zhf nj) 


( 12 ) 
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With regard to the shape of the functions ip only states of the same parity will give nonvanishing elements 
(r|s) so the states |0e0h), \0e2h), \le3h) etc. will be taken into account. To summarize in order to calculate 
the optical response of a wide Quantum Well it is necessery to solve the constitutive equation m using the 
matrix elements (r|s) and the potential matrix elements (HU. 


3 The parameters of the effective potential 


The further calculations require the estimation of parameters characterizing the effective potential m- 
We make the following assumptions: 1) The potential is isotropic, in analogy to the Coulomb potential in 
isotropic materials. The nanostructure anisotropy is included in the quasiparticles effective masses. This 
assumption leads to the equality u = v. 2) We assume the value S ~ 2R* {R* being the effective excitonic 
Rydberg energy for the given crystal); the exact value S will be established later. We determine the ground 
state energy of a hydrogen-like atom, where the interaction between the charges is given by (fTUl) . To this 
end we solve the Schrbdinger equation 


2 ^ \dr2 



Ip - SR*e-^^"-ip 


Eip. 


(14) 


Making use of the relation ^ = R*a*^ 


with the effective Bohr radius a*, we introduce scaled variables 


P = 


r 



vj = wa , 



transforming the eg. (1 1411 to Hip = eip with the Hamiltonian 


H = - 


\dp^ pdpj 


Se-^p\ 


(15) 


(16) 


The considered Schrbdinger equation will be solved by the variational method. Using the trial function 


iP = 


(17) 


we arrive at 


£:(A) = -X-S 
^ ’ 2 


\\ + wj 


3/2 


(18) 


By assuming the condition e = —1 and the vanishing derivative e' = 0, for any given value of A, one obtains 
a system of equations for two unknown quantities S and w. 



S 


A 


X + w 


3/2 


1 

s 



- 1 , 

0 , 


(19) 


and their values will be than used to determine the elements (fT^ . Looking for a solution which will reproduce 
the exact energy value £ = —1 we choose A = 0.34, S = 2.22, and w = 0.1. In order to compute the optical 
spectra we have to solve the system (|TT|) of coupled differential equations, but it will be easier to obtain the 
solutions by transforming the equations into linear algebraic equations. This can be done in the following 
way. Assume, for a moment, that the equation with indices (0,0) decouples from the remaining equations. 
Denoting Voooo = Vq we obtain the following equation for the amplitude Too 


Eg + Eoe + Eoh — hu — ir — -— 

2/X|| 


f d^ 1_^ 1 d^ \ 

\ dp2 p dp d(/)^ ) 



Yoo = Me 


Hp- Po) 

2 npo 


U(0|0). 


( 20 ) 


After rescaling the spatial variables in the effective excitonic Bohr radius the above equation becomes 
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Ap^ p dp p^ d(/)^ 
where now p denotes the scaled variable p/a*, and 


- aoe 


d/)o = ^(0e|0h), 

fi 27rpo 


^2 _ 
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Eg + Eoe + EqIi — fiU) — lE 
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Assuming the s-symmetry for the ground state, we first solve the Schrodinger equation 

d2 1 d 


dp^ p dp 


aoe 


Ip = £tp. 


( 21 ) 


( 22 ) 


(23) 


Using the variational method we solve above equation, using the trial function ip = e Denoting by H 

the left-hand-side operator, we compute the expression which should be minimized 


£(A)=<^ = A- 


Xvo 


{1P\1P) 


(A -I- tn)' 


The condition for the minimum yields 


where x = ^. The function 


(to — Vo)x^ -I- 2TOa; -|- to = 0, 


V'o(/0,(/) = 




(24) 


(25) 


(26) 


with the value of A obtained from the above equation can be considered as the normalized eigenfunction of 
the Schrodinger equation with the Hamiltonian 


d^ 1 d 

dp^ p dp 


Eo = —t~^ — :— — voe — Pp — ■cooooe 


(27) 


The index 0 denotes that this is the lowest energy state for the relative electron-hole motion with the assumed 
effective e-h interaction potential. Now we put doo(p) into the form 


Aoo(p) = Aipoip). 


(28) 


With regard to Hoipo = eoipo, where eo = e corresponds to the above estimated energy value, we obtain from 

(EU 


A = 


1 2p|| 

fcoo E £o Tt? 


MoE{Q\G)iPo{po). 


and the amplitude Y (p, Ze, Zh) has the form 


1 2aii 

Y (p, Ze, Zh) = 72 —- ^MoE{0e\0h)lpo (po) '0o(p)V'Oe (^e) V'o/i {zh) ■ 

^00 + ^0 a 

4 The solution of the constitutive equation 

Making use of the above calculated function tpo, we put the amplitude (H)) into the form 

N 

Y {P, Ze, Zh) = iJoip) X] i’ej{Ze)4’nh{Zh)Yjn, 

j,n=0 


(29) 


(30) 


( 31 ) 


where now Yjn are constant coefficients. Equation m takes now the form 
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\Ers + ^ ) i>o{p)Yrs - e tpo{p) VrsnjYnj = E-^^{er\hs)6 (p - Po) . (32) 

V ^^ 11 / 

After rescaling the spatial variable p —>■ p/a* we obtain from (l32ll the relation 

{krs + Pp) Mp)Yrs - e“^'’Vo(p) ^rsnjYnj = ‘^E^^{er\s)6 (p - Po), (33) 

which, using the quantities k/.^ = Vrsnj = can be written as 

(kls + Pp - ipo{p)Yrs + vooooe~'^P^tpo{p)Yrs 

-e~^p''ipo{p)'^VrsnjYnj = ^£1-^ (er |hs)(5 (p - po), (34) 

^ £ 27rpo 

and, in consequence, 

(fe^^ + eo) Yrs + Woooov ^—Yrs - ^VrsnjYnj = ^£1 A/q( er|hs)V'o (po) • (35) 

A -T VJ A “T tZ7 ^ fl 

nj 

We obtained a system of linear algebraic equations for the coefficients Ynj. Having them, we determine 
the amplitude Y (or amplitudes, when accounting the heavy- and light hole excitons H and L. Given the 
amplitude, we compute the polarization inside the quantum well and the electric field. For the further 
calculations we introduce dimensionless quantities Ws 

2Mo . . 

- = yrs-E 36 

eoehTTo* 

and arrived to the formula 

{krs Y Co) 3^rs + 'CqOOO T ' 3^rs T ' - V Vrsnjynj = (er|hs)V'o (po), (37) 

A -T ZD A ZD ' ^ ih 

nj 

where we used the relation 2^^ ~ with Alt being the transversal-longitudinal splitting energy, 

(see for example. [24)1. The described method can be used when we define the confinement energies fiuje, hujh 
and thus the parameters cth We will choose them to compare our theoretical results with the experimental 
findings of Miller et al.[I]. They obtained optical spectra for GaAs(Well)/Gao. 7 Alo. 3 As (Barrier) QWs of 
three thicknesses: £ = 51 ± 3,5 nm, £ = 32,5 ± 3.5 nm, £ = 33.6 ± 3.5 nm. It can be noticed the 
uncertainty in determining the well thickness. The confinement parameters are obtained as follows. We 
consider a symmetric QW with a rectangular confinement potential V 

V = £;g(Gao, 7 Alo, 3 As) - £;g(GaAs) = 482.8 meV, (38) 

see Tabled) The confinement potentials for electrons 14 and holes 14 are chosen as 

Ve = 0.85 V = 410.38 meV, 14 = 0.15 V = 72.42 meV. (39) 


Then we compute the lowest energy states in the QW with potentials 14,14- We follow the scheme from 
Ref. [29] where the lowest energies result from the equation 


ruB ) 


-f 1 


-1 

— cos^ P = 0, 


where the dimensionless parameters 6, 0q are defined as 


kL 1 [W L 2 mwVL^ l/HW£y 


(40) 
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and the index W means Well; the values a*, R* are appropriate for electrons and holes for the QW material, 
and are defined as 


R* 


4 

me 

2(47reoeb)2?i^ ’ 


* 


a 


h^jATreoeb) 

me^ 


(42) 


The below listed values are obtained when we insert in (1421) the appropriate effective masses: mg for 
and for and I^\\h,l are the in-plane reduced masses for the electron-hole pair and for 

the heavy- and light-hole exciton data. 


Table 1: Band parameter values for GaAs, AlAs, and Gao, 7 Alo, 3 As, AlAs data from [30], for Gao.rAlo.sAs 
by linear interpolation. Energies in meV, masses in free electron mass mo, 71,72 are Luttinger parameters 


Parameter GaAs AlAs Gao 7 AI 0 sAs 


me 

71 

72 

mh\\H 

mh\\L 

^^\\H 

M||l 

muzH 

mhzL 

E?* 

Hr 

Rl 

D* 

* 


e-b 


1519.2 

0.0665 

6.85 

2.1 

0.112 

0.210 

0.042 

0.05 

0.38 

0.09 

3.64 

4.3 

5.76 

15.78 

13.265 

9.97 

12.53 


3130 

0.124 

3.218 

0.628 

0.26 

0.386 


0.51 

0.22 

13.32 

19.35 

7.03 

4.84 

11.16 


2002 

0.084 


0.39 

0.13 


12.12 


First we determine the electron energy. For the further calculations we choose the well of GaAs thickness 
L = 51 nm. For the considered GaAs/Gao. 7 Alo. 3 As QW we have (Table |T]) mw = 0.0665 mo, ms = 
0.0840 mo. Using the values for GaAs from table |T] and substituting L = 51 nm into eq. dSD we obtain 
00 = 21.59. With this value we have 6 = 4.972 from (1401) and the lowest electron energy 

Eeo = 4i?: = 21.78 meV. (43) 

Quite analogous calculations can be performed for heavy- and light holes. For the heavy hole one obtains 

(44) 


a 1 jmhzHVh L 

t>0H = 7:1 


0'\\hR 


H “77 


Putting the appropriate data from TablejTjwe have Oqh = 20.61, 9h = 4.99 and the heavy-hole energy 


EhozH = EhoH = R*h = 4.23 meV. 

mhzH \ L J 

For the light hole Oql = 10.82, 6 ^ = 5.27 and the energy 

EhozL = EhoL = Rl = 17.20 meV. 

mhzL \ L J 


(45) 


( 46 ) 
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Thus the lowest confinement energy for the pair electron-heavy hole results 


EqzH — EeOz + -E'/tOzif — 21.78 -|- 4.23 — 26.01 meV 
and for the pair electron-light hole 

Eqzl = EeOz + EhozL = 21.78 + 17.20 = 38.98 meV. 

Now we identify the confinement energies with the lowest parabolic confinement energies: 

Tiujp^ „ tv-Oh 


Epn — 


Efio = 


and obtain the confinement parameters a 


OteO-H — 


meUJe 




I mhzH EhoH 


I rUe EeO 


= 3.08. 


= 3.07, 


(47) 

(48) 

(49) 

(50) 

(51) 


with analogous calculations for the light hole. For the pair electron-heavy hole (heavy-hole exciton) we 
obtain 


VboOO r,i^eHO‘hH)^lH 

UOOOOH = 


R 


H 


a 


eH 


■ a 


hH 


a 


eH 


a 


eH 


^hH 


^eH^hH 


-\- W 


a 


eH 


a 


hH 


- 1 - 1/2 


where 


aeH = afjae, 


(XhH = anO^hH- 


(52) 


(53) 


Making use of eqn. O, (EH), and putting w = 0.1, we obtain noooo = = 1-98. This value inserted into 

eq. (EH) gives X = 0.184 and A = 0.545, and from (1^ the lowest heavy-hole exciton energy cqh = —1.128. 
The lowest absorption peak observed in Ref. [1] corresponds to the energy 1535 meV, and the highest at 
about 1750 meV. Our calculations give the lowest heavy-hole exciton energy at 


Eg + if(eO) -|- E[h0) -|- £qhR*h — Eg E^q -|- E^qh + £ohRh ~ 1541 meV. 

The resonance at 1750 meV will be obtained for the state |e4h4), i.e. 

Eg + 2 ^ + ^4 -f —^ huJh + £ohRh — Eg + 9(Eeo + Ehon) + £ohRh ~ 1749 meV. 

The lowest resonance for the light-hole exciton is at energy 

Eg + E{e0) + i?(h0) -|- Eg = Eg + E^q + E^ql + SqlRl ~ 1553 meV. 
whereas for the state |e2h2) we have 


(54) 


(55) 


(56) 


1 


1 


Eg + ^2 -|- — J huje + ^2 -|- — j huih + £olRl — Eg -\- 5(ife0 + EhoL) + £olR*l ~ 1709 meV. (57) 

Thus we conclude that the resonances observed in Ref. [Tj come from the confinement states labeled by 
quantum numbers 0, 1, 2, 3, 4. 

As it follows from the relations (HH), (HH, and EH, the nonvanishing elements (er|hs) will be obtained 
for the confinement functions of the same parity, it means that either r = 2k, s = 2m; fc, m = 0,1,2 or 
r = 2k + 1, s = 2m + 1. The same holds for the potential matrix elements. With regard to this property we 














choose the following 13 electron-hole states with appropriate renumbering (both for heavy- and light-hole 
exciton) 


|e0h0) -i- |1), 

|elhl) ^ |2) 

|e2h2) ^ |3), 

|e3h3) ^ |4), 

|e4h4) -)> |5), 

|e0h2) ^ |6), 

|e0h4) |7), 

|elh3) |8), 

|e2h0) |9), 

|e2h4) |10), 

|e3hl) ^ 111), 

|e4h2) ^ |13), 

|e4h0) ^ |12), 


(58) 


where the notation means, for example 

\e2\)ff) = il;e 2 {ze)'>pho{zh), etc. 

The same operation is performed for energies for light and heavy hole excitons 


Eq -\- Eqr -\- Ehs — hui — lE 


Ejh, 


Eg EqT -|- Els ~ — ir —>■ Ejl, j — 1,2,..., 13, 


The potential matrix elements become now a square matrix 


Vrsnj = (erhs| exp[-'!; {ze - Zhf]enhj) Vji. 


(59) 

(60) 

(61) 


Using this notation we transform the equations (13711 into a system of linear equations for the 13 unknown 
quantities yj 

(fcy -|- eo) 3 ^j + uii-^—~ Vjnyn = bj 


‘ A -I- zu" 


A -I- tx7 


where 


bj = ^^(er|hs)V'o(po)- 


Equation (j62|l can be written in a matrix form 


where 3^, B are vectors 


Ay = B, 


y = {yi,...,yi3), 5 = (6i,...,5i3), 

and the matrix elements of A are defined as 

A 


(62) 

(63) 

(64) 

(65) 


+ ^0 
A 


A -I- Ci7 


A -I- t37 


(vii - Vjj), 

n + j- 


( 66 ) 


5 Results for GaAs/Gai_a;Ala;As parabolic Quantum Well and dis¬ 
cussion 

We have computed the optical functions of a GaAs/Gai_a;Ala;As parabolic Quantum Well with a chosen 
total thickness of 51 nm. The values of the relevant parameters are well known, and are given in Table 1. 
In our scheme the polarization inside the QW is related to the coherent amplitudes 


N 


Y{p,Ze,Zh) = ijjQ ^ |ejhn)Y,„, 

j,n=0 


(67) 
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by the relation 


N 

P{z) = 2MotpQ{po) y] \ejhn){z), (68) 

j,n=0 

with the notation 

|ejhn)(z) = ipej{z)il’hniz). (69) 

Having the polarization, we compute the mean dielectric susceptibility 

N 

X = 7rebV'o(po)y]3^«(lK)A/2 (70) 

e=o 

A/2 

where (1|^)a/2 = x / K)(C)dC, A = Having the susceptibility, one can compute, using the appropri- 

-A/2 

ate boundary conditions, the optical functions (reflectivity, transmission, and absorption). We choose the 
absorption, which is related to the effective dielectric function by the formula 


a = —Imi/eT+X (71) 

being the dielectric constant of the QW material. Now we can compare the theoretical absorption spectra 
obtained by (1711) with the luminescence spectra from Ref.[T]. We have computed the absorption coefficient 
for the described above Wide parabolic QW of the thickness 51 nm. The first step was to determine the 
coefficients S, w satisfying the equations CH). Then, by using the potential partition (1391) and the formerly 
obtained value u, we have computed the potential matrix elenmts Vrsnj and the matrix elements (rls). 
Assuming a certain value of the coherence radius po, we have determined the lowest excitonic eigenfunction 
//’o- Finally, taking a certain value of the damping parameter T, we have solved the constitutive equation m, 
obtaining the coherent amplitudes. From the amplitudes we have computed the mean dielectric susceptibility 
dZOl) and the absorption coefficient (ED. The results for the real and imaginary part of the mean susceptibility 
of the considered QW are displayed in Fig. [D The parameters used in the calculations are listed in the figure 
caption. The arrows indicate the positions of absorption maxima from Ref. [!]• The good agreement of 
theory and experiment (both in positions of maxima and their oscillator strengths) can be seen. In general, 
we observe 17 resonance peaks, from which 15 can be identified with those observed in experiment. The 
detailed comparison with peaks enumerated by rising energy is shown in the Table [2] We have chosen the 
parameters to obtain the best fit to the experimental results of Ref. [1] . The accuracy of the optimal choice 
of the effective potential parameters and damping can be tested in the following way. We have computed 
the total fitting error for the first 13 maxima as a function of the parameters S, zu, F and v. The results 
are shown on the Fig. [2] (a) and (b). We learned that the positions of the absorption maxima is mainly 
affected by the values w and S. One can see that the change the values of these parameters stretches the 
whole spectrum, causing a linear shift of the peak position, as shown on the Fig. (c). When using the 
value S = 2.6, we obtain w Ri 0.154, which represents a local minimum of fitting error. The assumed 
value of u = 0.5 is also a good choice. For the global minimum at u = 1.2, some parts of the absorption 
spectrum became negative, which was deemed unphysical. As expected, small values of F have no effect on 
the location of the peaks. For significant values of F, some peaks become indistinguishable, which is seen as 
a sudden jump in the fitting error. The selected parameter values gave the theoretical maxima close to the 
experimental values with mean error of less than 3.5 meV and enabled to identify the electron-hole states. 

In the next step we tried to fit the experimental line shapes (oscillator strengths). We have observed 
that variations of the coherence radius change substantially the lineshapes. The best fit was obtained for 
Pol = 0.17 a*j^,pQH = 0.1 a^. It can be also verified that the increase of the damping parameter T results 
the lowering of the oscillator strength. 
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Energy (meV) 


Figure 1: The real and imaginary part of the mean QW susceptibility for the heavy-hole (H) and light-hole 
(L) exciton. The parameters used in calculations are zu = 0.154, S = 2.6, v = 0.5, F = 0.5 and the coherence 
radii poL = 0.17 a\,poH = 0.1 a*jj, respectively. The electron-hole states and their energies are assessed 
(indexed by 1,..., 11) and the corresponding maxima from [I] are indicated below, with the heights of the 
bars indicating the oscillator strengths. 


Table 2: The identification of the electron-hole states 


Number of maximum 

State description 

Nearest maximum from Ref. [1]. 

1 

Elh (1533 meV) 

Elh (1530 meV), Ell (1536 meV) 

2 

Ell (1546 meV) 

E13h (1550 meV) 

3 

E6h (1566 meV) 

E2h (1565 meV), E21 (1567 meV) 

4 

E7h (1583 meV) 

E24h (1583 meV) 

5 

E2h (1588 meV) 

E24h (1583 meVh E3h (1595 meV) 

6 

E8h (1618 meV) 

E35h (1613 meV) 

7 

E21 (1632 meV) 

E4h (1626 meV), E46h (1640 meV) 

8 

E9h (1636 meV) 

E46h (1640 meV) 

9 

E61 (1636 meV) 

E46h (1640 meV) 

10 

E3h (1642 meV) 

E46h (1640 meV) 

11 

E91 (1654 meV) 

E5h (1655 meV) 

12 

ElOh (1671 meV) 

E57h (1671 meV) 

13 

E4h (1688 meV) 

E6h (1686 meV) 

14 

Ellh (1696 meV) 

E6h (1686 meV), E68h (1700 meV) 

15 

E71 (1704 meV) 

E68h (1700 meV) 
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Energy [meV] total fitting error [meV] total fitting error [meV] 



(a) Total fitting error 





(b) Total fitting error as a function of v and T. 




(c) The effect of the parameters S and zu on the position of the hrst three heavy hole exciton peaks. 


Figure 2: The choice of the optimal calculation parameters. 
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6 Conclusions 


We have developed a simple mathematical procedure to calculate the optical functions of wide parabolic 
quantum wells. Our procedure describes the optical properties of a QW, taking into account the Coulomb 
interaction between electrons and holes. Our treatment includes anisotropic properties of the QW, and takes 
into account coherence of the electron-hole pair with the radiation field. The presented method has been used 
to investigate the optical functions of GaAs/Gai_a;Al 2 ;As parabolic Quantum Well for the case of radiation 
incidence parallel to the growth direction and it shows an excellent agreement with the experimental data, 
explaining the number and the positions of the absorption maxima. The justification of the choice of effective 
potential parameters and the damping constant is also presented. 
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